Investigating the Cosmic-Ray Ionization Rate in the Galactic 
Diffuse Interstellar Medium through Observations of H3 

Nick Indriolo^'^, Benjamin J. McCall^'^ 
ABSTRACT 

Observations of H.^ in the Galactic diffuse interstellar medium (ISM) have led 
to various surprising results, including the conclusion that the cosmic-ray ioniza- 
tion rate (C2) is about 1 order of magnitude larger than previously thought. The 
present survey expands the sample of diffuse cloud sight lines with observa- 
tions to 50, with detections in 21 of those. Ionization rates inferred from these 
observations are in the range (1.7 ± 1.3) x lO'^^ s'^ < C2 < (10.6 ± 8.2) x 10"^^ 
s^^ with a mean value of C2 = (3-5l;3;o) x 10^^^ s^^ Upper limits {3a) derived 
from non-detections of B.^ are as low as (2 < 0.4 x 10~^® s~^. These low upper- 
limits, in combination with the wide range of inferred cosmic-ray ionization rates, 
indicate variations in (2 between different diffuse cloud sight lines. A study of 
(2 versus iVn (total hydrogen column density) shows that the two parameters 
are not correlated for diffuse molecular cloud sight lines, but that the ionization 
rate decreases when A^h increases to values typical of dense molecular clouds. 
Both the difference in ionization rates between diffuse and dense clouds and the 
variation of (2 among diffuse cloud sight lines are likely the result of particle prop- 
agation effects. The lower ionization rate in dense clouds is due to the inability 
of low-energy (few MeV) protons to penetrate such regions, while the ionization 
rate in diffuse clouds is controlled by the proximity of the observed cloud to a 
site of particle acceleration. 
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INTRODUCTION 



The rich chemistry responsible for producing many of the small molecular species ob- 
served in the interstellar medium (ISM) is thought to be driven primarily by ion-molecule 
reactions ( IWatsonI Il973l : iHerbst fc Klempererl Il973l ). To initiate this reaction network, a 
source of ionization is required that operates in both diffuse molecular (n ~ 100 cm~^, 
T ~ 70 K) and dense molecular (n ~ 10^ cm~^, T ~ 30 K) clouds. As absorption by atomic 
hydrogen severely attenuates the flux of photons capable of ionizing many of the most abun- 
dant species in the molecular ISM (e.g., H2, He, O), cosmic rays are the dominant ionization 
mechanism in such regions. This makes the cosmic-ray ionization rate a vital parameter in 
modeling interstellar chemistry. 

Several estimates of the cosmic -ray ionization rate have been made over the past 50 years 
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Following its initial detection in the ISM (iGeballe fc Okalll996l). Hj" has wid elv been 
regarded as the most direct tracer of the cosmic-ray ionization rat e (lDalgarndl2006[) due to a 
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been used to constrain the ionization rate in each of these environments. Here we present the 
most comprehensive survey of Hjj" in d iffuse molecular clou d sight lines to date, expanding 
and improving upon our analysis from llndriolo et al.l (120071 ). 



2. H+ CHEMISTRY 



As mentioned above, the chemistry associated with H^ in the ISM is rather simple. Hj] 
is formed in a two-step process, beginning with the ionization of H2 by cosmic rays. 



H2 + CR ^ H+ + + CR' 



(1) 
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and quickly followed by a reaction of with H2, 

H+ + H2 H+ + H. (2) 
Some amount of will be destroyed by dissociative recombination with electrons, 

H+ + e- ^ H + H, (3) 
or charge transfer to atomic hydrogen, 

H+ + H ^ H2 + H+, (4) 

before it can form H^, but these processes are slow compared to reaction ([2]). Overall, 
reaction ([T]) is the rate-limiting step (it is many orders of magnitude slower than reaction 
(j2])) and can be taken as the formation rate of Hj]". Photoionization of H2 also occurs, but 
is negligible compared to ionization by cosmic rays. This is because the ultraviolet and soft 
X-ray photons capable of ionizing H2 {E > 15.4 eV) are severely attenuated in the outer 
layers of diffuse clouds by atomic hydrogen, which has a lower ionization potential (13.6 eV). 
While higher energy photons (e.g., hard X-rays and 7-rays) can penetrate diffuse molecular 
clouds, the lower ionization cross sections at these energies and smaller flux of such photons 
make their contribution to the ionization rate minimal. As a result, Hjj' should be formed 
primarily through the ionization of H2 by cosmic rays. 

The primary destruction mechanisms for Hjj' are dependent on the environment under 
consideration. In diffuse molecular clouds, Hj]" is predominantly destroyed via dissociative 
recombination with electrons, 

H+ + e" ^ H2 + H or H + H + H. (5) 

In dense clouds however, where the electron fraction is much lower, is destroyed by proton 
transfer to neutrals such as CO, O, and N2: 

H+ + CO ^ H2 + HCO+, (6) 

H+ + CO ^ H2 + HOC+, (7) 

H+ + O ^ H2 + 0H+, (8) 

H+ + N2 ^ H2 + HN+. (9) 

These ion-molecule reactions demonstrate how Hj]" essentially drives the chemistry in the 
ISM by generating molecular ions which can then go on to form even more complex species. 
Because the chemistry surrounding Hjj' is so simple, and because its formation arises directly 
as a result of ionization by cosmic rays, observations of H;^ can be a powerful probe of the 
cosmic-ray ionization rate. 
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METHODS 



3.1. Observations 



In diffuse molecular clouds, conditions are such that only the 2 lowest energy levels of 
— the (J, fC) = (l,l) para and (1,0) ortho states — are expected to be significantly popu- 
lated. As a result, observations that probe transitions arising from these 2 states will trace 
the entire content of B.^ along a line of sight. We targeted 3 such transitions — the -R(l, 1)", 
R{1, 0), and R{1, 1)' — which are observable at 3.668083 fim, 3.668516 /im, and 3.715479 /im, 
respectively. As the weak Hj]" lines are located near various atmospheric CH4 and HDO 
lines, precise removal of the telluric lines is vital for the purpose of accurately determining 
column densities, and is facilitated by the observation of telluric standard stars. Because 
the i?(l, 1)" and R{1, 0) transitions are easily observed simultaneously (being separated by 
about 4 A), and because together they probe the entire population of H3 , they are the most 
frequently observed transitions. 

Although observations of H.^ have become more commonplace, as evidenced by the pub- 
lication list in Section [H t hey are still rather difficult. Given the transition dipole moments 
(see, e.g.. lGoto et ahlbood . and references therein) and a typical relative abundance of about 
10~^-10~^ with respect to H2, H3 absorption lines in diffuse molecular clouds are usually 
about 1-2% deep, and thus require a signal-to- noise ratio (S/N) greater than about 300 on 
the continuum to achieve 3a detections. As a result, background sources must have L-band 
(3.5 /im) magnitudes brighter than about L = 7.5 mag to be feasible targets for 8 m (or 
L = 5.5 mag for 4 m) class telescopes. Ideally, these background sources should also be 
early-type stars (OB) such that the continuum is free of stellar absorption features. A list 
of all background sources targeted in our study, along with select line-of-sight properties, is 
presented in Table [H Additionally, to adequately sample the intrinsically narrow interstellar 
absorption lines requires high spectroscopic resolving power {R > 20,000). 

Observations used in this study were taken bet ween 1998 and 2010 with a variety of 
instrur nent / telescope comb inations, including CGS4 (jMountain et al.lll990r) at UKIRT, NIR- 



SPEC flMcLean et al. 



19981) at Keck. P hoenix flHinkle et al.ll2003r i at both KPNO and Gemini 



South, and CRIRES (IKaufl et al.ll2004j ) at VLT. A log of all Hj^ observations in diffuse molec- 
ular clouds utilized herein — including those previously pubhshed — can be found in Table 



2.4 of llndriold (120 111 ). For specifics regarding instrument confi gurations, observing meth- 



ods a n d the lik e , the r eader is referred to th e following sour c es: iMcCall et al.l (119981 . 
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3.2. Data Reduction 



Beginning with raw data frames, our reduction procedure employed standard IRAF0 
routines commonly used in spectroscopic data reduction. Upon extracting one-dimensional 
spectra, d ata were tranf erred to Igor Prci whe re we have macros w ritten to complete the 
reduction ( IMcCalll l200ll ). Spectra published in llndriolo et al.l ( 120071 ) were reprocessed here 
using a method that better accounts for bad pixels in the CGS4 array. A full description 
of the new reduction procedure used for CGS4 data — and all of our reduction procedures 
applying to B.^ data in general — is presented in llndriold ( 1201 ll ). 



4. UPDATED RESULTS 



Fully reduced spectra are presented in Figures [THIS] • Of the 42 sight lines presented, 
15 show absorption lines, w hile 27 do not. Combining these results with the spectra 
presented in iMcCall et al.l (120021 ). our entire survey covers 50 diffuse molecular cloud sight 



lines, with Hg detections in 21 of those. 



As mentioned above, all of the spectra presented in llndriolo et al.l ( 120071 ) were re- 
processed for this study. In most cases the differences between the resulting spectra are 
minimal, but spectra presented herein, along with any results derived from them, should 
be regarded as superseding any previously published results. There are 4 sight lines in 
particular where further discussion of new versus old spe ctra is warranted: BD -14 5037, 
HD 154368, ( Per, and X Per. In llndriolo et al.l (120071 ) we reported a detection of Hjj" 
toward BD -14 5037 — albeit marginal at best — but with the updated reduction procedure 
this sight line is now cons idered to be a non-detection (see Figure |6] herein versus Figure 3 
from llndriolo et al.l (120071) ). Spectra toward HD 154368, ( Per, and X Per were reported in 



Crabtree et al.l (l201ll ). but we later realized that an important step in the reduction process 
was omitted for these 3 sight lines; spectra taken at the A and B nod positions had been 
combined based solely on pixel number without accounting for any potential shift in the 
dispersion direction. For both HD 154368 and C Per the effects of re-processing the data 
were mini mal, and newly deriv ed column densities and para fractions agree with those re- 
ported in ICrabtree et al.l (120111 ) within uncertainties. For X Per though, individual spectra 
had relatively low S/N and a large number of frames were shifted in t he dispersion di r ection . 



thus the new reduction is significantly different from that reported in ICrabtree et al.l (120111 ) . 



^http:/ /iraf.noao.edu/ 
^http://www.wavenietrics.coni/ 
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By us ing the new X Per results (see Table [2]) in the analysis presented by ICrabtree et al. 



( 120 111 ), the conclusions in that study are actually strengthened, as the erroneous X Per data 
used therein consistently made that sight line the only outlier in the reported trends. 



5. ANALYSIS 

5.1. Extraction of Column Densities and Upper Limits 

Absorption features due to were fit with Gaussian functions in order to determine 
equivalent widths, velocity FWHM, and interstellar gas velocities. The fitting procedure 
uses the functional form of a Gaussian where the area (as opposed to amplitude) is a free 
parameter, and includes a fit to the continuum level; i.e.. 



Iq — exp 

W\/Tf 



A- Ar 



w 



(10) 



and was developed in Igor Pro. The free parameters here are the continuum level, Iq, the area 
of the Gaussian, A, the central wavelength of the Gaussian, Aq, and the line width, w (n.b., 
the width used here is related to the "standard" Gaussian width, a, by w'^ = 2cr^). As the 
continuum level has already been normalized, an area determined using this fit is by definition 
an equivalent width, Wx. In the case of the i?(l, 1)" and i?(l,0) lines, both absorption 
features are fit simultaneously and a single best-fit continuum level is found. Uncertainties 
on the equivalent widths (SWx) and continuum level {61) — both at the la level — are 
output by the fitting process. To estimate the systematic uncertainties due to continuum 
placement, the continuum level was forced to Jq + SI and Jq — 61 and the absorption lines re- 
fit. Variations in the equivalent widths due to this shift are small compared to those reported 
by the fitting procedure and so have been ignored (i.e., (j{Wx) = SWx). Assumir ig optically 



thin a bsorption lines and taking transition dipole moments and wavelengths from lGoto et al. 



( I2OO2I . and references therein) , column densities are derived from equivalent widths using the 
standard relation: 

where N{J, K) is the column density in the state from which the transition arises, h is 
Planck's constant, c is the speed of light, A is the transition wavelength, and is the 
square of the transition dipole moment. 

In cases where H3 absorption lines are not detected, upper limits to the equivalent width 
are computed as 

Wx < (y\^^M~, (12) 
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where a is the standard deviation on the continuum level near the expected H3 lines, Apix is 
the wavelength coverage per pixel, and ATpix is the number of pixels expected in an absorption 
feature. Upper limits to the column density are then determined via equation flTT]) . Column 
densities, equivalent widths, velocity FWHM, and interstellar gas velocities for all of the 
sight lines presented in Figures [T]-[T3] are reported in Table [2l 



5.2. Determination of the Cosmic-Ray Ionization Rate 

As discussed in Section [21 the chemistry associated with Hj]" in diffuse molecular clouds 
is rather simple. Assuming steady-state conditions, the formation and destruction rates of 
can be equated. In the simplified version of diffuse cloud chemistry — where destruction 
occurs only via dissociative recombination with electrons — the result is 

C2n(H2) = A;en(H+)ne, (13) 

where C2 is the ionization rate of H2, n's are number densities, and fee is the Hj]"-electron 
recombination rate coefficient (see Table [3]). Substituting the electron fraction (defined as 
Xe = Ue/n-fi, where nn = n(H) + 2n(B.2)) into equation (fT3!l and solving for the ionization 
rate gives 

C2 = Kx.n^^r^. 14 
n(H2) 

Observations cannot measure changes in these parameters along a line-of-sight, so we assume 
a uniform cloud with path length L and constant Xe, fee, ^h, and n(Hjj")/n(H2). In this 
case, the Hjj' and H2 densities can by definition be replaced with iV(Hjj')/L and N{^2)jL^ 
respectively, such that 

gives the cosmic-ray ionization rate in a diffuse molecular cloud. 

While equation (IT^ only considers one formation and destruction mechanism for H^, 
comparison to a more complete chemical reaction network shows that it is a robust approx- 
imation given diffuse cloud conditions. Assuming steady-state for where destruction by 
electron recombination and charge transfer to protons (reactions [3] and H]) are considered, 
and steady-state for Hjj" where destruction by proton transfer to CO and O (reactions O El 
and [HD is accounted for, gives the equations 

C2n(H2) = ^en{^t) + Agp(H)n(H+) + A|2p(H2)n(H+), (16) 

A[2p(H2)n(H+) = Kn^ni^X) + A;co«(CO)n(H+) + fcon(0)n(H+). (17) 
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Solving for the ionization rate and making similar substitutions as before results in 



C2 



-nn [keXe + /ccoa;(CO) + A;oa;(0)] 



where relevant rate coefficients are given in Table [3], /cco 
is the fraction of hydrogen nuclei in molecular form. 



2^ 



1 



A[g]+ /|7], and = 2n(H2)/nH 



The cosmic-ray ionization rate as a function of electron fraction is plotted in Figure 
Hlfor both equations (ITSll and (|T8l) . In all four panels the thick solid line shows the linear 
relationship given by equation f|T5l) . while the various other curves show the ionization rate 
determined using equation f|T8|) and different values of a;(CO), a;(0), and /hj- In the first 
panel a;(CO) is set to 10-^ 10-^ and 10"^ with a;(0) = 10"* and /h^ = 1. In the second 
panel x(0) is set to 10"^ 10-^ 10"^ and 3 x 10"^ with x(CO) = 10"^ and /h^ = 1. In the 
third panel /h^ is set to 0.67, 0.5, and 0.1 with x(CO) = x{0) = 10~«. 



The first and second panels show that equation (fT8|) only differs significantly from equa- 
tion (fT5|) when Xg < 10~^, and for high fractional abundances of CO and O (> 10~^). These 
deviations occur when proton transfer to CO and O come to dominate over electron recom- 
bination as the primary destruction pathways for H3 . The third panel shows that affects 
the ionization rate for all values of Xe, but that rather low molecular hydrogen fractions are 
necessary to significantly alter the inferred value of ^2- With /hj = 0.5 (i.e., twice as many 
H atoms as H2 molecules) the value output by equation ( |T8|) is about 1.6 times that from 
equation f[T^ . while for /hj = 0.67 (i.e., equal number of H and H2) the value output by 
equation (fT8|) is about 1.3 times that from equation (fTSl) . This deviation is caused by the 
larger relative abundance of atomic hydrogen destroying before it can form Hjj", and is 
represented by the final term in equation (fTSll . The destruction of by electron recombi- 
nation does not play a major role in the chemical network used here, and its infiuence can 
only be seen in the slight deviation between equations (ITBl) and (ITSl) when Xe ~ 10~^. 

Although some of the sight lines studied herein ha ye molecular hydrogen fractions of 



fn-, = 0.2 a s determined from observations of H and H2 (ISavage et al.l 119771 : iRachford et al. 



2OO2I . l2009l ). it must be remembered that these are line-of-sight fractions, while equation 
(ITS]) requires local fractions. As it is expected that there will be atomic gas along a line of 
sight that is not associated with a cloud containing H2, the hne-of-sight value of always 
underestimates the value of within a diffuse molecular cloud. It is generally assumed that 
the interior of a diffuse molecular cloud has conditions where 0.67 < < 1 (i-e. somewhere 
between half and all of the hydrogen is in molecular form), such that reaction (jlj) will not 
be very important. 



Usi ng average d i ffuse molecular cloud c onditions (x(0) = 3 x 10 ^ ICartledge et al. 



(2004 


); 


Jensen et al. 


( 


2005) 



10-6 ISonnentrucker et al.l fl2007f ): /h^ > 0.67) in 
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equation (fT8|) produces the dashed curve in the las t panel of Figure [T^ and at Xe = x{C^) = 
1.5 X 10~^ ( jCardelli et al.lll996t ISofia et al.l 120041 ) the resulting ionization rate is only 1.33 
times larger than that inferred from equation f[T5]l . This demonstrates that equation (llSp is 
a good approximation to equation f lTS]) in the regions we study, and so we adopt equation 
( IT5|) in computing ^2- 

With equation ( !T5|) . the cosmic-ray ionization rate (C2) can be determined from the 
column density (A^(H3 )), H2 column density (A^(H2)), total hydrogen density (nn), electron 
fraction (xe), and H^-electron recombination rate coefficient (ke)- The electron fraction can 
be approximated by the fractional abundance of assuming that nearly all electrons are the 
result of singly photoionized carbon. In sight lines where C"*" has been observed, the measured 
value of x{C'^) is used for Xe- For all other sight lines, t he average fractiona l abundance 
measured in various diffuse clouds, x(C^) ~ 1.5 x 10~^ ( ICardelli et al.l Il996l : ISofia et al. 



20041 ). is adopted for Xg. Uncertainties in Xe are assumed to be ±20%, i.e., ±3 x 10 ^. The 



Hj]'-electron recor nbination rate coefficient has now been meas ured in multiple laboratory 



experiments (e.g.. iMcCall et al.ll2004j : iKreckel et al.ll2005l . 120101 ) with consistent results and 
is presented in Table El When available the spin temperature of H2, Tqi, and its uncertainty 
are used in calculating ke] otherwise, an average value of 70 ± 10 K is adopted. The total 
hydrogen number density is difficult to determine, but vari ous studies have estimated tt-h 
using a rotational excitation analysis of observed C 2 lines (ISonnentrucker et al.l 120071 ). an 
analysis of H and th e J = 4 level of (|Juralll975l ). or a thermal pressure analysis of fine 
structure lines of C I ( iJenkins et al.lll983l ). Number densities from these studies are presented 
in Table H] when available. For sight lines without estimated densities, the rough average 
value of riH = 200 cm~^ is adopted. In all cases, uncertainties in nn are assumed to be ±50% 
of the reported values. Absorption lines from electronic transitions of H2 have been observed 
in the UV along m any of the sight lines in this study (from which the aforementioned values 
of Tqi are derived; ISavage et al.lll977l : iRachford et al.ll2002l . |2009|). In sight lines where H2 
has not been observed, two other methods were used to estimate A^(H2). The preferred 
method uses column densities of CH determined from observa t ions i n combination with the 
relation A^(CH)/iV(H2) = 3.5^11 x 10~^ from Isheffer et al.l feoosh . In sight lines where 
neither H2 nor CH has been observed, A^(H2) is estimated from th e color excess , E(B — V), 
using the relation A^h ~ E{B — V)5.8 x 10^^ cm~^ mag~^ from iBohlin et al.l ( 119781 ). and 
assuming /hj ~ 2A^(II2)/A'h = 0.67. Molecular hydrogen column densities determined both 
from observations and estimates are presented in Table HI 

While Tabled gives individual column densities for the lowest lying ortho and para levels 
of H3 , values of A^(H3 ) in Table IH are equal to the sum of the column densities in the (1, 1) 
and (1,0) states (i.e., A^(Hjj") = A^(l, 1) ± A^(1,0)). In cases where observations of both 
the i?(l, 1)" and i?(l, 1)' lines produce independent values of A^(l,l), a variance weighted 
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average, 

n n 

N = Y.^N,/a^)/Y,{llal), (19) 

i=l i=l 

is used to determine A^(l, 1). Because there are only 2 independent measurements of A^(l, 1), 
the uncertainty of the weighted average, 5'[A^(1, 1)], is determined using an unbiased estima- 
tor of a weighted population variance for small samples. This is given by 

V " 



where 



Wi = \, Vi = ^ Wi, V2 = ^ Wi, 
i=i i=i 



and /i* is the weighted average determined from equation f|T9|) . The uncertainty in the 
total H3 column density is then computed as usual by adding a[N{l, 1)] and cr[A^(l,0)] in 
quadrature. Upper limits to the H3 column density should be taken as 3o"[A^(H3 )]. 

Using equation (fT5|) . and taking the values described above and in Table HJ cosmic-ray 
ionization rates are inferred for all diffuse molecular clouds where Hj]" observations have been 
made. These values of C2 are presented in column 10 of Table HI with uncertainties in column 
11. As for the Hj|" column densities, upper limits to the cosmic-ray ionization rate should be 
taken as 3o"(C2)- 

For one sight line however, that toward NGC 2024 IRS 1, a different analysis is used 
because recent observations suggest that the interstellar material is more likely dense than 
diffuse (T. Snow 2011, private communication). In this case, values appropriate for dense 
clouds (xe = 10"'', = 1) are adopted, effectively simplifying equation (ITSl) to 

C2 = ^^^H [fccox(CO) + A;ox(0)] . (21) 

The density and temperature are also set to average dense cloud values (nn = 10^ cm~^, 
T = 30 K). The large CO column density (iV(CO) = 1.26 x 10^^ cm'^; T. Snow 2011, private 
communication) results in x{CO) = 1.28 x 10~^, demonstrating that most of the carbon is in 
molecular form, and thus validating the low value o f Xe assumed a bove. Oxygen abundances 



are typically about two times carbon abundances ( lLoddersll2003l ). such that half of all O is 
expected to be in the form of CO. We assume the remainder to be in atomic form, and use 
a;(0) = a;(CO) in equation ( |2T]) . The resulting upper limit on the cosmic-ray ionization rate 
toward NGC 2024 IRS 1 is C2 < 4.2 x 10^^'' s^^ This lower value is expected in dense cloud 
conditions, as will be discussed below, and is excluded from our analysis of the distribution 
of ionization rates in diffuse clouds. 
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5.3. The Distribution of Cosmic- Ray Ionization Rates in Diffuse Clouds 



With ionization rates and upper hmits inferred for 50 diffuse cloud sight lines, it is 
interesting to study the distribution of our entire sample. Cosmic-ray ionization rates with 
la uncertainties and 3a upper limits are plotted in Figure [T5l In order to find the probability 
density function for the cosmic-ray ionization rate in our sample, all of these data points are 
combined. Each sight line with an Hjj" detection is treated as a gaussian function with mean 
and standard deviation equal to the ionization rate and uncertainty reported in columns 10 
and 11 of Table m Each sight line with a non-detection of Hj]" is treated as a gaussian function 
with mean equal to zero and standard deviation equal to the uncertainty in column 11 of 
Table ID The cosmic-ray ionization rate cannot be negative, and it is reasonable to assume 
that in diffuse clouds there is some "floor" value below which G does not fall caused by the 
continuous diffusion of particles throughout the Galactic disk (IWebberlll998l ). We set this 
floor value to 10~^^ s~^ — a rou gh average of ionization rates computed from proposed local 



interstellar proton spectra (e.g.. ISpitzer fc Tomaskdll968l : IWebber 



19981 ) — and truncate all 



Each function is then re-normalized so that the area 
under the curve above the floor ionization rate is equal to 1. The truncated gaussians are 
added together, and the flnal distribution is again normalized to have an area equal to 1. 

The resulting probability density function, /(C2), is plotted in the top panel of Figure 
[TBI Because of the logarithmic x-axis, it is difficult to tell "by eye" which ionization rates 
are most probable. However, this is easily seen by scaling /(C2) by (2, as shown in the 
middle panel of Figure [T6l By integrating (2!' {(2), we flnd the mean value of the cosmic-ray 
ionization rate to be 3.5 x 10~^^ s~^, marked by the vertical dashed line. The minimum 
range of ionization rates (in log space) that contains 68.3% (gaussian la equivalent) of the 
area under /(C2) is bounded by 5.0 x 10~^^ s~^ on the left and 8.8 x 10~^^ s~^ on the right, 
and is shown by the shaded region. The bottom panel of Figure [16] shows the cumulative 
distribution function, F(C2), which gives the probability that the ionization rate is below a 
particular value. 



■^Setting the floor ionization rate to zero does not significantly change our results. 
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CORRELATIONS BETWEEN LINE-OF-SIGHT PROPERTIES 



6.1. The Column Density 

Given that the formation process starts with H2, one might expect that abundances 
of the 2 species should be correlated. Figure [T71 shows A^(Hj^) versus A^(H2) for the entire 
sample of sight lines studied he rein. Also included a re column densities in si ght lines near 



the supernova remnant IC 443 f|lndriolo et al.ll2010r), as we ll as in dense clouds (IMcCall et al. 



I999I : iKulesal I2OO2I : iBrittain et al.l 120041 : iGibb et al.ll2010f ). It is clear that for diffuse cloud 
sight lines where H3 has been detected, relative abundances tend to cluster about 10~^. A 
re- arrangement of equation (fT5|) . 



C2 



A(H2) KXeUYl' 

suggests that scatter about a central value of ^(Hg )/A(H2) is hkely due to variations in ^2 
and electron density (xenn) between sight lines. 

The much lower value of ^(Hg )/A(H2) in dense clouds (10^*^-10^^) is the result of a 
higher density, different destruction partner, and lower ionization rate. Replacing electrons 
with CO as the dominant destruction partner leaves N{B.^)/N{B.2) inversely proportional to 
kcox{CO)nii. The rate coefficient for destruction of by CO is about 2 orders of magnitude 
slower than by electrons (see Table BI), the density of dens e clouds is about 2-3 orders of 
magnitude larger than diffuse clouds (ISnow fc McCallll2006l ). and the fractional abundance 
of electrons in diffuse clouds should be about equal to that of CO in dense clouds. This 
means that (2 must be about 0-2 orders of magnitude lower in dense clouds than in diffuse 
clouds to produce the obser ved results, and several studies of ionized species i n dense clouds 



supp ort such a trend (e.g., Ivan der Tak fc van Dishoeckl I2OOOI : iKulesal |2002[ iHezareh et al. 
2008h . 



6.2. The Cosmic-Ray Ionization Rate and Location 

Beyond studying the cosmic-ray ionization rate in and of itself, it is insightful to search 
for correlations between (2 and various other sight line parameters. One of the most fun- 
damental relationships to study is that between ^2 and location within the Galaxy. The 
cosmic-ray ionization rate is plotted against Galactic longitude in the bottom panel of Fig- 
ure [m The top panel is a map in Galactic coordinates and the middle panel gives the 
distance to background sources such that the reader can easily trace the ionization rates 
in the bottom panel to an actual on-sky position. There are no apparent gradients in Fig- 
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ure [IHl and comparisons of C2 with both hehocentric distance and Galactocentric radius of 
background sources do not show any trends, thus suggesting that the mechanism responsible 
for variations in the cosmic-ray ionization rate does not operate on large scales. To study 
variations of (2 on small spatial scales, we investigate some of the closest sight line pairings 
from our data. 



6.2.1. HD 168607 & HD 168625 

The sight lines with the smallest angular separation in our survey are HD 168607 and 
HD 168625, separated by only 1^1. The Hj]" spectra for both targets are very similar (see 
Figure [2]), suggesting that the absorption arises from a common cloud that must be in front 
of the closer target (HD 168607) at 1100 pc. At this distance, I'.l corresponds to an on-sky 
separation of 0.35 pc. Ionization rates for these two sight lines are consistent with each other 
within uncertainties, demonstrating uniformity of C2 on the smallest spatial scales probed by 
our observations. 



6.2.2. Cyg 0B2 Association 



Three of our target sight lines use stars in the Cyg 0B2 association as background 
sources: Cyg 0B2 5, Cyg 0B2 12 , and Cyg 0B2 8A. Taking the OB association to be at 
1.7 kpc (jTorres-Dodgen et al.lll99l[ ). the on-sky distance between No. 5 and No. 12 is 2.6 pc, 
between No. 5 and No. 8A is 4.9 pc, and between No. 12 and No. 8A is 3.8 pc. Spectra 
toward No. 5 and No. 12 are somewhat similar, with comp onents at about 6 km s~^ and 
12 km s~^ LSR in a few molecular species (IMcCall et al.ll2002l ). While the ionization rates in 
these two sight lines differ by about a factor of 3, the inability to meaningfully separate the Hjj' 
absorption features into the different velocity components makes this analysis very uncertain, 
and for Cyg 0B2 5 and Cyg 0B2 12 have overlapping la uncertainties. Cyg 0B2 8A, 
although not far removed from the other two sight lines, do es not show Hj" a bsorption (see 
Figure [T2|) . and may only contain the 6 km s~^ component (jSnow et al.ll2010l ). Still, the 3(T 
upper limit on (,2 toward No. 8 A is consistent with the ionization rates inferred for the other 
two sight lines, and nothing conclusive can be drawn from this region. 

A complementary method of tracing the cosmic-ray flux relies on gamma-ray emission 
from 7r° decay, and observations of the Cygnus X region by Fermi LAT are sug gestive of 
a large flux of recently accelerated hadronic cosmic rays (lAckermann et al.ll201ll ) near the 
Cyg 0B2 association. This is one of a few select regions where both gamma-rays and H^ 
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have been observed, and presents the opportunity for a rather unique study based on both 
cosmic ray tracers. However, this analysis is beyond the scope of the current paper, and 
such a study would best be served by an expanded survey of Ht that probes a larg er extent 
of the gas traced by the gamma-ray emission presented in lAckermann et al.l (l201l[ ) . 



6.2.3. Per 0B2 Association 
The targets C, Per, X Per, and a Per are all in the Perse us 0B2 association ( although 



the latter two may only be in the association by projection; Ide Zeeuw et al.lll999[ ). X Per 
and C Per show very similar spectra (see Figure [T]) with Hj]" absorption lines at 7-8 km s~^ 
LSR, and the ionization rate inferred for each sight line is nearly identical. At a distance of 
301 pc, the angular separation of 0?88 corresponds to 4.6 pc in the plane of the sky. Among 
our targeted sight lines o Per is the next closest to these two, and, assuming a distance of 
283 pc, is 10.5 pc away from C, Per and 13.2 pc away from X Per. H3 is not detected toward 
Per, and the 3(T upper limit inferred for C2 in this sight line is about equal to the lower la 
bounds on the ionization rates toward X Per and C, Per. While Per is about 20 pc closer 
than C Per, the similar E{B — V) values for both sight lines make it unlikely that we do not 
see HjI" simply because it is behind the target. Although not definitive, these results strongly 
suggest variations in the cosmic-ray ionization rate on a size scale of order 10 pc. 



6.2.4. Ophiuchus-Scorpius Region 

Many of our target sight lines probe the nearby Ophiuchus-Scorpius region, and some 
of the most closely spaced include Sco, p Oph D, HD 147889, and x Oph. The first three 
are all separated by 2-3 pc (assuming a distance of 136 pc), with x Oph about 12-14 pc 
away. This is one of a few regions where is not detected in any of the target sight lines, 
and aside from p Oph D where the low continuum S/N level (~ 175) of the spectrum limits 
our analysis, the 3a upper limits on (2 are in the range 3-19 x 10~^^ s~^. These are some 
of the smaller inferred upper limits on (2, and likely pertain to the material closest to the 
sun probed in our survey as background sources are only 100-200 pc away. It is interesting 
then, that these upper limits are consistent with the ion ization rate of a few times 10^^^ s~^ 



predicted by the local interstellar cosmic-ray spectrum (jWebberlll998l ). 
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6.3. The Cosmic-Ray Ionization Rate and Hydrogen Column Density 



Given that there is no evidence for large-scale gradients in the cosmic-ray ionization rate, 
it is possible that properties of the observed clouds themselves are responsible for variations 
in and will show some correlation with the cosmic-ray ionization rate. A parameter with 



which (2 has been predicted to vary is the hydrogen column density, A^h fjPadovani et al. 
20091 ). This is because the energy spectrum of cosmic-rays is expected to change with 
depth into a cloud. Lower-energy particles — those most efficient at ionization — will lose 
all of their energy to ionization interactions in the outer regions of a cloud, leaving only 
higher-energy particles to ionize the cloud interior. As such, the ionization rate in a cloud 
interior should be lower than in a cloud exterior. Because lower- energy particles can operat e 
through a larger portion of clouds with lower column densities ( jCravens fc Dalgarndll978l ). 
it is expected that the inferred (2 will decrease with increasing A^h- A plot of (2 versus A^h 
(and equivalent E{B — V)) is shown in Figure [T9l In sight lines with multiple distinct velocity 
components, the total line-of-sight A^h has been divided by the number of components to 
better approximate the size of each indivi dual cloud. In cluded in Figure ^\ are data from 4 

(2002), 5 dense cloud sight lines observed 



it lines observed in by 



Ku 



esa 



dense cloud sig 

in H^^CO"'" bv Ivan der Tak fc van Dishoeckl (j2000l ). and 1 dense cloud sight line observed 
in HCNH"*" by Hezareh et al. ( 20081 ) for the purpose of extending the relationship to much 
higher A^h- There does not appear to be a correlation between (2 and A^h when only the 
diffuse molecular cloud sight lines are considered. However, when the dense cloud ionization 
rates (most of which are of order a few times 10^^'' s^^) are included, it seems that (2 does 
decrease in sight lines with higher column densities. 



7. DISCUSSION 

Given that (2 does not decrease with increased A'h among diffuse cloud sight lines, but 
does when switching from the diffuse to dense cloud regime, the following conclusions can be 
drawn. First, the cosmic rays that are primarily responsible for ionization in diffuse clouds 
must be able to propagate entirely through such clouds. This could easily be don e by 10 



MeV protons, which have a range of Rn ^ 2 x 10^^ cm~^ (jCravens &: Dalgarndll978l ). about 
equal to the largest values of A'h in diffuse clouds studied here. However, a cloud with a 
large line-of-sight column density does not necessarily have a large column density in the 
plane of the sky. Even lower-energy particles then (e.g., 2 MeV protons with Rn ^ 10^^ 
cm~^) could potentially cause ionization through the entire extent of an observed cloud. 



Second, the cosmic rays that are responsible for most of the ionization in diffuse clouds 
must not be able to penetrate to the interiors of dense clouds. Again, something like 2- 
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10 MeV protons fit this picture well. Most of the sight lines where Hjj' is observed in dense 
clouds use embedded objects as background sources. Regions surrounding these sources have 
been mapped in various molecular species (via emission at radio wavelengths), implying large 
column densities in the plane of the sky in addition to the large column densities observed 
along the line of sight. This increases the likelihood that much of the material being probed 
is deep within the observed cloud. 

Taken together, the higher ionization rates inferred in diffuse clouds and lower ionization 
rates inferred in dense clouds suggest that differences in (2 can be explained by the inabil- 
ity of low-energy cosmic rays to penetrate large columns of material. However, u nder the 
assum ption that the cosmic-ray spectrum is uniform throughout the Galactic disk (jWebber 
19981 ). this does not adequately explain differences in the cosmic-ray ionization rate amongst 
diffuse clouds. Instead, it would seem that the cosmic-ray spectrum — at least for particles 
in the energy range most efficient at ionization — must vary in space. 

Despite the previously held assumption of a uniform cosmic-ray spectrum, it should not 
be surprising that the fiux of low-energy particles varies across the Galaxy. Given a hydrogen 
density of 1 cm~^, a 2 MeV proton will only travel about 320 pc (not necessarily in a straight 
line) before losing all of its energy, meaning that any point that is a few hundred parsecs 
away from a site of particle acceleration will not experience the same fiux of 2 MeV protons 
as a point that is much closer to an acceleration site. For low-energy particles to even enter 
a diffuse cloud then, the cloud must be relatively close to a site of particle acceleration. 

Regardless of where most low-energy cosmic rays are accelerated, it is possible that 
differences in (2 among diffuse cloud sight lines can be attributed to the distance between 
a cloud and acceleration site. Sight lines that probe material in close proximity to an 
acceleration site should show high ionization rates, while those that probe material farther 
away should show lower ionization rates. Unfortunately, defining acceleration sites and 
actually computing physical distances between those and the sample of observed clouds is 
difficult at best. 

There are a few special cases, however, where the ionization rate can be inferred in 
molecular clouds known to be in close pro ximity to a site of cosmic-ray acceleratio n such as 



a supernova remnan t. Shocked gas (e.g. , iHuang et al. 



(1720 MHz) masers (IClaussen et al.lll997l : lHewitt et al. 



■19861 : iDickman et al.lll992f ) and OH 
20061 ) observed toward the supernova 



remnant IC 443 demonstrate that the expanding shell is p hysically interactirig with a fore 



ground cloud, and ganima-ray observatiori s of the remnant (lAlbert et al.ll2007l : lAcciari et al. 



2OO9I : lAbdo et al.ll2010l : iTavani et al.ll2010t ) show a signature indicative of 7r° decay that sug- 
gests the cloud is experiencing a large fiux of energetic protons. In llndriolo et al.l (120101 ) 
we presented observations of Hjj' in 6 sight lines that pass through molecular material near 
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the supernova remnant IC 443. For the 2 sight hnes where Hg was detected we inferred 
ionization rates of C2 = 2.6lj^'9 ^ 10^^^ s^^ and C2 = l-6^?'2 x 10~^^ s^^ — much higher than 
the mean value reported above — while in the remaining sight lines Sa upper limits on (,2 
were consistent with typical diffuse cloud values. Not only do these findings support the 
hypothesis that proximity to a site of particle acceleration controls the cosmic-ray ionization 
rate, but they also demonstrate variations in C2 on length scales of a few pc. 



8. SUMMARY 

The survey of Hjj" in diffuse molecular clouds now covers 50 sight lines, with detections 
in 21 of those. Cosmic-ray ionization rates (and upper limits) are inferred in all of these sight 
lines. Where Hjj" is detected, ionization rates are in the range (1.7 ± 1.3) x 10~^^ s~^ < C2 < 
(10.6 ± 8.2) X 10^^^ s~^. Accounting for upper limits, the mean value of the ionization rate 
is C2 = (3.5I3 0) X 10^^^ s^^. This is about 1 order of magnitude larger than ionization rates 
inferred from previous observations of other molecular species (e.g., OH and HD), although 
recent observations of 0H+ and H20^ with Her schel are suggestive of high ionization rates 



as well (0.6 x 10"^^ s^^ < Ch < 2.4 x 10"^^ s'^ iNeufeld et al.ll2010[ ). The lowest 3a upper 



limits found for sight lines where is not detected are about (2 < 0.4 x 10~^^ s~^. Together, 
the wide range of inferred ionization rates and the low upper limits allude to variations in 
(2 between different diffuse cloud sight lines. 

Comparisons of (2 with various line-of-sight properties (including Galactic latitude. 
Galactic longitude and heliocentric distance) show no clear trends. This suggests that vari- 
ations in the cosmic-ray ionization rate between different sight lines are caused not by large 
scale, but local effects. A comparison of C2 with the total hydrogen column density (Ah) 
shows no strong correlation when only diffuse clouds are considered. When ionization rates 
inferred for dense cloud sight lines with much higher Ah are included though, (2 is seen to 
decrease with increased Ah. This correlation is expected as low-energy cosmic rays — those 
most efficient at ionization — will lose all of their energy before reaching the interiors of 
dense clouds. The lack of a correlation amongst diffuse cloud sight lines suggests that the 
particles primarily responsible for ionization must be able to completely penetrate diffuse 
clouds, while the difference in ^2 between diffuse and dense clouds indicates that these same 
particles must be stopped in the outer layers of dense clouds. 

Given that total column density — at least in the case of diffuse clouds — is unable 
to explain variations in the ionization rate, the concept of a cosmic-ray spectrum that is 
uniform throughout the Galaxy must be reconsidered. Instead, the flux of low-energy (MeV) 
cosmic rays is most likely controlled by proximity to local sites of particle acceleration. Our 
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recent observations showing; a high ionization rate toward the supernova remnant IC 443 



( llndriolo et al.ll2010[ ) support this conjecture. 
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Table 1:: Science Targets 

E{B - V) Distance 



Object 


HD number 


(mag.) 


Rcf. 


(pc) 


Rcf. 


K Cas 


2905 


0.33 


1 


1370 


11 


HD 20041 


20041 


0.72 


2 


1400 


2 


HD 21483 


21483 


0.56 


2 


440 


2 


HD 21389 


21389 


0.57 


2 


940 


2 


HD 21856 


21856 


0.19 


3 


502 


11 


40 Per 


22951 


0.27 


1 


283 


11 


Per 


23180 


0.31 


2 


280 


2 


BD +31 643 


281159 


0.85 


2 


240 


2 


C Per 


24398 


0.31 


2 


301 


2 


X Per 


24534 


0.59 


2 


590 


2 


e Per 


24760 


0.09 


3 


196 


11 


^ Per 


24912 


0.33 


2 


470 


2 


62 Tau 


27778 


0.37 


2 


223 


2 


HD 29647 


29647 


1.00 


2 


177 


2 


a Cam 


30614 


0.30 


2 


820 


2 


NGC 2024 IRS 1 




1.69 


4 


450 


4 


Ori 


41117 


0.45 


2 


1000 


2 


HD 47129 


47129 


0.36 


1 


1514 


12 


HD 53367 


53367 


0.74 


2 


780 


2 


HD 73882 


73882 


0.70 


5 


1000 


13 


HD 110432 


110432 


0.51 


5 


300 


14 


Sco 


147084 


0.74 


4 


270 


11 


p Oph D 


147888 


0.47 


2 


136 


2 


HD 147889 


147889 


1.07 


2 


136 


2 


X Oph 


148184 


0.52 


1 


161 


5 


IJ. Nor 


149038 


0.38 


3 


1122 


12 


HD 149404 


149404 


0.68 


1 


417 


11 


COph 


149757 


0.32 


2 


140 


2 


HD 152236 


152236 


0.68 


5 


600 


15 


HD 154368 


154368 


0.78 


5 


909 


16 


WR 104 




2.30 


6 


2600 


17 


HD 168607 


168607 


1.61 


7 


1100 


7 


HD 168625 


168625 


1.86 


8 


2800 


18 


BD -14 5037 




1.57 


4 


1700 


4 


HD 169454 


169454 


1.12 


2 


930 


2 


W 40 IRS lA 




2.90 


9 


400 


9 


WR 118 




4.57 


6 


3900 


6 


WR 121 




1.75 


6 


2600 


6 


HD 183143 


183143 


1.28 


7 


1000 


7 


P Cygni 


193237 


0.63 


7 


1800 


7 


HD 193322A 


193322A 


0.40 


3 


830 


19 


HD 229059 


229059 


1.71 


2 


1000 


2 
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Table 1:: (continued) 



Object 


HD number 


E{B-V) 
(mag.) 


Ref. 


Distance 
(pc) 


Ref. 


HD 194279 


194279 


1.22 


7 


1100 


7 


Cyg 0B2 5 




1.99 


7 


1700 


7 


Cyg 0B2 12 




3.35 


7 


1700 


7 


Cyg 0B2 8A 




1.60 


10 


1700 


7 


HD 204827 


204827 


1.11 


2 


600 


2 


HD 206267 


206267 


0.53 


2 


1000 


2 


19 Cep 


209975 


0.36 


1 


1300 


20 


A Cep 


210839 


0.57 


2 


505 


2 


1 Cas 


218376 


0.25 


2 


339 


2 



Re ferences: (1) llVIcCall et all (Eolol): (2) iThorburn et al.l (l2003h : 



Rachford et al 



(3) ISavaee et al.l (119771) : (A) Ivan Dishoeck fc Blackl (119891): (51 



(|2009l) : C6l lConti fc Vaccal (Il990l): (7) 



ShuDing et al.. (1999): (10) |S 
J2007 ):d2') IPruscione et al 



now et al. (201 



Chevalier fc Ilovaiskv 



McCall et al. 



(I2OO2L and references th erein) ; (8) iNota et all (1 1996h : 



11) 



van Leeuwen 



1994); (13) 
1998|); (15) 



Nichols fc Slavin 



Ritchev et al 



2004); 



2Q11D; 



(14)^ 

(16) iMaiVApellaniz et a P (l2004h: (17) iTuthill et al.l (I2OO8I): (18) 



Pasauah et al.l|2002l ): (19) [Roberts et al.l|2O10l ): ("20) lsheffer et al 
(|200 " 
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Table 2:: Absorption Line Parameters 







fLSR 


FWHM 






N{J, K) 


a{N) 


Object 


Transition 


(km 


(km s^^) 


(10-^ fim) 


(10-® /J.m) 


(10^^ cm-2) 


(10" cm-2) 


HD 20041 


i?(l, 1)" 


2.5 


11.0 


2.49 


0.45 


10.3 


1.87 




i?(l,0) 


-0.9 


10.1 


3.03 


0.41 


7.67 


1.04 


HD 21389 


R{1, 1)" 


0.2 


13.3 


1.20 


0.19 


4.98 


0.79 




il(l,0) 


0.8 


12.9 


1.46 


0.19 


3.69 


0.47 


C Per 


R(l, 1)" 


8.0 


10.3 


0.99 


0.11 


4.09 


0.46 




R(1,0) 


8.4 


8.7 


0.86 


0.10 


2.17 


0.25 




i?(3,3)' 




10 




0.10 




0.37 


X Per 


_R(1, 1)" 


7.0 


11.7 


1.16 


0.20 


4.81 


0.81 




i?(l,0) 


6.3 


9.4 


1.00 


0.17 


2.53 


0.43 


HD 29647 


R(l, 1)" 


8.0 


8.6 


1.62 


0.21 


6.72 


0.88 




i?(l,0) 


8.2 


11.0 


2.06 


0.26 


5.21 


0.67 


HD 73882 


_R(1, 1)" 


5.9 


3.9 


1.44 


0.21 


5.97 


0.86 




K(1,0) 


5.7 


3.2 


1.16 


0.19 


2.94 


0.48 




R{i,iy 


5.4 


3.5 


1.34 


0.15 


6.15 


0.69 


HD 110432 


R{1, 1)" 


-3.8 


6.9 


0.74 


0.06 


3.08 


0.24 




ii(l,0) 


-3.3 


7.5 


0.83 


0.07 


2.11 


0.17 






-3.1 


8.1 


0.69 


0.06 


3.15 


0.28 


HD 154368 


R{1, 1)" 


5.4 


5.8 


1.57 


0.28 


6.51 


1.16 




R(1,0) 


4.6 


5.1 


1.13 


0.25 


2.86 


0.63 


HD 168607 


R{1, 1)" 


20.9 


7.9 


1.07 


0.13 


4.44 


0.53 




i?(l,0) 


22.4 


8.9 


0.85 


0.14 


2.16 


0.34 


HD 168625 


R{1, 1)" 


19.5 


9.5 


1.58 


0.19 


6.55 


0.78 




il(l,0) 


22.0 


9.1 


1.59 


0.19 


4.02 


0.48 


HD 169454 


R(l,l)" 


4.5 


10.6 


0.80 


0.12 


3.32 


0.51 




R(1,0) 


4.0 


11.7 


1.01 


0.12 


2.56 


0.30 




R(l,l)' 


3.8 


8.7 


0.77 


0.18 


3.53 


0.84 


W40 IRS lA 


_R(1,1)" 


5.9 


8.9 


4.19 


0.64 


17.4 


2.66 




i?(l,0) 


7.4 


7.8 


4.48 


0.60 


11.3 


1.51 




ii(l,l)' 


6.1 


4.5 


2.80 


0.71 


12.8 


3.26 


HD 229059 


R(l, 1)" 


3.8 


9.3 


5.10 


0.25 


21.2 


1.02 




i?(l,0) 


3.5 


9.7 


6.42 


0.26 


16.2 


0.65 


HD 204827 


_R(1,1)" 


-0.8 


11.0 


3.11 


0.53 


12.9 


2.18 




ij(l,0) 


-0.9 


8.7 


2.41 


0.51 


6.10 


1.29 


A Cep 


_R(1,1)" 


-2.9 


10.7 


1.04 


0.23 


4.31 


0.95 




i?(l,0) 


-1.4 


13.5 


1.29 


0.27 


3.26 


0.67 


K Gas 


R(l, 1)" 




15 




0.19 




0.77 




R(1,0) 




15 




0.19 




0.47 


HD 21483 


_R(1,1)" 




15 




0.47 




1.97 




R{1,0) 




15 




0.47 




1.20 


HD 21856 


R(l,l)" 




15 




0.46 




1.92 




R{1,0) 




15 




0.46 




1.17 


40 Per 


R{1, 1)" 




15 




0.18 




0.74 




i?{l,0) 




15 




0.18 




0.45 


o Per 


R(l, 1)" 




10 




0.10 




0.41 




R(1,0) 




10 




0.10 




0.25 


BD +31 643 


_R(1,1)" 




15 




0.95 




3.95 




i?(l,0) 




15 




0.95 




2.41 


e Per 


R{1, 1)" 




10 




0.08 




0.33 




il(l,0) 




10 




0.08 




0.20 


? Per 


R(l,l)" 




10 




0.11 




0.44 




R(1,0) 




10 




0.11 




0.27 


62 Tau 


_R{1, 1)" 




15 




0.37 




1.53 




i?(l,0) 




15 




0.37 




0.93 


a Cam 






15 




0.30 




1.26 
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Table 2:: (continued) 



Object 


""LSR 

Transition (km s^-"^) 


FWHM 

(km s^-"-) 


(10-^ tim) {IQ-^ fim) 


NiJ, K) 
(10" cm-2) 


a{N) 
(10" cm~2) 




/?(1,0) 


15 


0.30 




0.77 


NGC 2024 IRSl 


iJ(l, 1)" 


10 


0.26 




1.09 




R(1,0) 


10 


0.26 




0.66 


HD 47129 


i?(l,l)" 


15 


0.60 




2.48 




i?(l,0) 


15 


0.60 




1.51 


HD 53367 


fl(l,l)" 


5 


0.12 




0.50 




i?(l,0) 


5 


0.12 




0.30 




R(l,l)' 


5 


0.07 




0.32 


o Sco 


fl(l,l)" 


10 


0.08 




0.35 




i?(l,0) 


10 


0.08 




0.21 


HD 147888 


i?.(l,l)" 


10 


0.46 




1.92 




i?(l,0) 


10 


0.46 




1.17 


HD 147889 


i?.(l,l)" 


10 


0.18 




0.74 




i?(l,0) 


10 


0.18 




0.45 


X Oph 


iJ(l,l)" 


5 


0.11 




0.45 




R(1,0) 


5 


0.11 




0.27 


fi Nor 


ii;(l,l)" 


5 


0.12 




0.49 




R(1,0) 


5 


0.12 




0.30 


HD 149404 


i?(l,l)" 


5 


0.23 




0.94 




R(1,0) 


5 


0.23 




0.57 


C Oph 


iJ(l,l)" 


10 


0.08 




0.31 




R(1,0) 


10 


0.08 




0.19 


HD 152236 


i?.(l,l)" 


5 


0.14 




0.59 




fl(l,0) 


5 


0.14 




0.36 


BD -14 5037 


fl(l,l)" 


10 


0.23 




0.95 




i?(l,0) 


10 


0.23 




0.58 


HD 193322A 


ii;(i,i)" 


15 


0.37 




1.55 




i?(l,0) 


15 


0.37 




0.95 


Cyg OB2 8A 


ii:(i,i)" 


15 


0.98 




4.07 




i?(l,0) 


15 


0.98 




2.48 


HD 206267 


ij(i,i)" 


10 


0.30 




1.25 




R(1,0) 


10 


0.30 




0.76 


19 Cep 


^(1,1)" 


15 


0.26 




1.10 




i?(l,0) 


15 


0.26 




0.67 


1 Gas 


ii:(i,i)" 


15 


0.29 




1.21 




R(1,0) 


15 


0.29 




0.74 



Notes: Columns 6 {a(W\)) and 8 {ct{N)) arc Icr uncertainties on the equivalent width and column density, respec- 
tively. The FWHM values in the case of non-detections are appropriate for the resolution of the instrument used in 
each observation, and are used in computing the value of A/'pix required in equation I I12I I. 



Table 3:: Rate Coefficients for Reactions Involved in Ht Chemistry 



Reaction 




Rate Coefficient (cm'^ s ^) 


Reference 


H+ 


+ H2 - 


^H+ + H 


2.08 X 10-9 


1 




+ e- - 


^ H + H 


A[3]= 1.6 X 10-*^(r/300)-0-^3 


2 


H+ 




H2 + H+ 


/g]= 6.4 X 10-1° 


3 




+ e- - 


■> products 


f^=ke = -1.3 X 10-8 + 1.27 X 10-6r-0-48 


4 




+ C0 


H2 + HCO+ 


A[g]= 1.36 X 10-9(T/300)-0-i42gxp(3.41/r) 


5 




+ C0 


H2 + HOC+ 


^= 8.49 X 10-10 (r/300)°-066i exp(-5.21/r) 


5 




+ 0^ 


H2 + 0H+ 


Ajg]= ko = 1.14 X 10-9(T/300)-°-i56exp(-1.41/r) 


5 




+ N2 - 


-> H2 + HN+ 


%]= 1.8 X 10-9 


6 



Notes: Numerical subscripts on the rate coefficients refer to the corresponding reaction 
numbers given in the text. The Hjj'-electron recombination rate coefficient, is referred 
to as ke throughout the text. The rate coefficient for destruction via proton transfer to CO 
used in the text, kco, is equal to /i[g]+ The rate coefficient for destruction via proton 
transfer to O, 1^ is referred to as ko throughout the text. All temperature dependent 
coefficients are in terms of the gas kinetic temperature, T, except for k^, which is in terms 
of the electron temperature, Tg. However, it is expected that electrons and H2 should be 
thermalized to the gas kinetic temperature via collisions, so we make no distinction between 
Tg , Toi , and T in our calculations. 

Re ferences: (1] iTheard fc H untressI ^197^)■. ^2 ) iMitchelll kmd): ^3) iKarnas et al.l (|l97d ): 
(4) iMcCall et all (|2004l ): fs;) iKlip penstein et al.l toid ): r6) lR,akshitl (Il982l ) 



Table 4:: Cosmic- Ray Ionization Rates 



Object 


A'(HJ) 
(10" cm-2) 


<7[Af(H+)] 
(10" cm-2) 


Reference 


Af(H2) 
(10^0 cra-'^) 


o-[Ar(H2)] 
(10^" cm"2) 


H2 
Reference 


(cm-3) 


Density 
Reference 


C2 

{W-^^ s^i) 


(10^^*^ s 


HD 20041 


23.7 


5.13 


1,2 


11.4 


6.24 


5 


200 




9.48 


7.59 


HD 21389 


8.67 


0.92 


1 


5.57 


3.06 


5 


200 




7.11 


5.54 


C Per 


6.26 


0.52 


1 


4.75 


0.95 


6 


215 


10 


5.55 


3.18 


X Per 


7.34 


0.92 


1 


8.38 


0.89 


7 


325 


10 


5.85 


3.54 


HD 29647 


11.9 


1.11 


1 


22.9 


12.5 


5 


200 




2.38 


1.85 


HD 73882 


9.02 


0.50 


3 


12.9 


2.39 


7 


520 


10 


9.71 


5.57 


HD 110432 


5.22 


0.17 


3 


4.37 


0.29 


7 


140 


10 


3.86 


2.10 


HD 154368 


9.37 


1.32 


1 


14.4 


3.99 


7 


240 


10 


4.19 


2.62 


WR 104 


23.2 


1.54 


2 


44.7 


22.3 


8 


200 




2.37 


1.76 


HD 168607 


6.60 


0.63 


1 


17.4 


8.47 


5 


200 




1.73 


1.28 


HD 168625 


10.6 


0.92 


1 


16.3 


7.95 


5 


200 




2.96 


2.17 


HD 169454 


5.93 


0.34 


1 


16.6 


8.37 


5 


300 


10 


2.45 


1.83 


W40 IRS lA 


26.9 


3.54 


1 


56.3 


28.2 


8 


300 


11 


3.27 


2.45 


WR 118 (5) 


29.4 


7.07 


2 


43.3 


21.7 


8 


200 




3.10 


2.41 


WR 118 (47) 


30.8 


5.73 


2 


45.5 


22.7 


8 


200 




3.10 


2.36 


WR 121 


22.0 


2.60 


2 


34.0 


17.0 


8 


200 




2.96 


2.21 


HD 183143 (7) 


11.3 


2.98 


2 


4.86 


2.38 


5 


200 




10.6 


8.23 


HD 183143 (24) 


13.5 


2.81 


2 


7.91 


3.85 


5 


200 




7.82 


5.93 


HD 229059 


37.4 


1.20 


1 


65.7 


51.2 


5 


200 




2.60 


2.47 


Cyg OB2 5 


24.0 


3.29 


2 


15.2 


7.39 


5 


225 


10 


8.13 


6.03 


Cyg OB2 12 


34.3 


5.89 


2,4 


80.0 


69.1 


5 


300 


10 


2.93 


3.04 


HD 204827 


19.0 


2.54 


1 


20.9 


10.2 


5 


450 


10 


9.32 


6.92 


A Cep 


7.58 


1.17 


1 


6.88 


0.48 


7 


115 


10 


2.84 


1.61 


K Cas 




0.91 


1 


1.88 


0.20 


6 


200 






1.78 


HD 21483 




2.30 


1 


11.9 


5.91 


5 


200 






0.88 


HD 21856 




2.25 


1 


1.10 


0.12 


6 


200 






8.42 


40 Per 




0.87 


1 


2.92 


0.53 


6 


80 


12 




0.57 


Per 




0.47 


1 


4.09 


0.79 


6 


265 


10 




0.85 


BD +31 643 




4.63 


1 


12.4 


4.39 


7 


200 






1.67 


e Per 




0.38 


1 


0.33 


0.07 


6 


15 


13 




0.36 


5 Per 




0.51 


1 


3.42 


0.53 


6 


300 


13 




1.98 


62 Tau 




1.79 


1 


6.23 


1.12 


7 


280 


10 




2.06 


a Cam 




1.48 


1 


2.17 


0.26 


6 


150 


13 




2.17 


Ori 




2.75 


2 


4.90 


1.05 


9 


200 






2.79 


HD 47129 




2.91 


1 


3.51 


0.58 


6 


200 






3.81 


HD 53367 




0.44 


1 


11.1 


1.44 


9 


200 






0.20 


o See 




0.41 


1 


19.1 


9.81 


5 


225 


10 




0.11 


HD 147888 




2.25 


1 


2.97 


0.61 


9 


215 


10 




9.25 


HD 147889 




0.87 


1 


28.8 


14.0 


5 


525 


10 




0.36 


X Oph 




0.52 


1 


4.26 


1.02 


6 


185 


10 




0.64 


IX Nor 




0.58 


1 


2.77 


0.49 


6 


200 






1.09 



Table 4:: (continued) 



Object 
HD 149404 

C Oph 
HD 152236 
BD -14 5037 

P Cygni 
HD 193322A 
HD 194279 
Cyg OB2 8A 
HD 206267 
19 Cep 
1 Gas 
NGC 2024 IRS 1 



(10^ 



a[Nin+)] 
(10^3 cm-2) 


T-r-i- 

H3 

Reference 


Af(H2) 
(10^0 cm-2) 


a[N(H2)] 
(lO^O cm-2) 


H2 
Reference 


nil 
(cm"'') 


Density ^2 
Reference (lO^^'^ s"!) 


^^^2) 
(lO^i** s 


1.10 


1 


6.17 


0.53 


9 


200 




0.87 


0.37 


1,2 


4.49 


0.42 


6 


215 


10 


0.39 


0.69 


1 


5.34 


1.32 


9 


200 




0.63 


1.11 


1 


38.6 


23.6 


5 


200 




0.13 


1.84 


2 


3.14 


1.61 


5 


200 




2.67 


1.82 


1 


1.21 


0.18 


6 


200 




6.52 


4.59 


2 


23.7 


11.9 


8 


200 




0.88 


4.77 


1 


9.26 


4.58 


5 


200 




2.35 


1.47 


1 


7.18 


0.59 


7 


200 




0.97 


1.28 


1 


1.21 


0.18 


6 


200 




4.60 


1.41 


1 


1.40 


0.19 


6 


200 




4.18 


1.27 


1 


49.0 


24.5 


8 


10000 




0.14 



Notes: The sight lines toward WR 118 and HD 183143 have 2 distinct velocity components observed in Ht bv lMcCall et al.l ( 2002 ). 
and each of these components is labeled with a number in parentheses that corresponds to its average LSR velocity. When no reference 
for TT-H was available the value was set to 200 cm~^, an average density for diffuse molecular clouds. Upper limits to the cosmic-ray 
ionization rate should be taken as 3 times a(C2)- 

H+ References: (1) this work; (2) iMcCah et all |2002l ): (3) ICrabtree et"al] |201ll ): (4) iMcCah et all |l998h 

H 2 R efere nces: (5) ca l culate d fro m CH column d e nsitie s using the relation iV(CH)/7V(H2) = 3.5^1,1 x 10"* from ISheffer et al 



I2OO8I ): (6) ISavage et al 
E{B-V)5.8x 10^1 — 



(I1977I) 



mag" 



; (7) iRachford et alj (120021) : (8) calculated from E{B - V) (see Table [I| using the rela t ion JV u 
-1 fromlioh lin et all (llQTSi). and assuming 2i V(H2) /A''H = /ho = ;67: o T lRachford et al.| (|2009( ) 
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Fig. 1. — Spectra targeting the R{1, 1)" and -R(l, 0) lines of Hj]" in the sight hnes toward HD 
20041, HD 21389, C Per, X Per, and HD 29647. Vertical hnes mark the expected positions of 
absorption lines given gas velocities along each line of sight. Gaps in spectra are where the 
removal of atmospheric features was particularly poor. Observations were made with CGS4 
at UKIRT. 
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Fig. 2. — Same as Figure [T] except for the sight hnes toward HD 168607, HD 168625, HD 
204827, and HD 229059. 
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Fig. 3. — Same as Figure [1] except for the sight hnes toward HD 169454 and W 40 IRS lA. 



- 33 - 




Fig. 4. — Same as Figure [3] except targeting the R{1, 1)' transition. 
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Fig. 5. — Same as Figure [T] except for the sight hnes toward o Per, e Per, ^ Per, NGC 2024 
IRS 1, and HD 47129. 
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Fig. 6. — Same as Figure [T] except for the sight hnes toward o Sco, HD 147888, HD 147889, 
C Oph, BD -14 5037, and HD 206267. 
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Fig. 7. — Same as Figure [T] except for the sight hnes toward HD 53367, HD 73882, and HD 
110432. Observations were made with CRIRES at VLT. 
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Fig. 8. — Same as Figure [7] except targeting the R{1, 1)' transition. 
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Fig. 9. — Same as Figure [T] except for the sight hnes toward x Oph, Nor, HD 149404, HD 
152236, and HD 154368. Observations were made with Phoenix at Gemini South. 
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Fig. 10. — Same as Figure [T] except for the sight hnes toward HD 21483, 40 Per, and A Cep. 
Spectra combined observations using CGS4 at UKIRT and NIRSPEC at Keck. 
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Fig. 11. — Same as Figured] except for the sight hnes toward HD 21856, BD +31 643, 62 Tau, 
and a Cam. Observations were made with NIRSPEC at Keck. 
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Fig. 12.— Same as Figure [U] except for the sight hnes toward HD 193322A, Cyg 0B2 8A, 
19 Cep, 1 Cas, and k Cas. 
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Fig. 13. — Spectrum covering the i?(3,3)' line of toward C Per taken with CGS4 at 
UKIRT. The vertical line marks the expected position of the absorption line. Absorption 
from the (3,3) metastable state is not expected in the diffuse ISM, and thus far has only been 
detec ted in the Galactic Center JCoto et al.ll2002l . l2008l . l201ll : IOka et al.lboOsl i lGeballe fc Oka 

2ninl ). 
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Fig. 14. — Cosmic-ray ionization rate (C2) as a function of the electron fraction (xe) for simplified chemistry (equation [15]) and more 
complete chemistry feauationlTS)). In all cases na = 200 cm" -3, 7V(H+)/Af(H2) = 10"^ and T = 70 K. The thick sohd line shows the linear 
approximation given by equation (jTS)). The first panel shows the result of varying a;(CO), with a;(0) = 10^® and — 1; the second panel 
shows the result of varying a;(0), with a::(CO) = 10^* and = 1; the third panel shows the result of varying /h2, with x{CO) ~ 10^* 
and x(0) = 10"^; the last panel shows equation ([H]) for average diffuse ISM conditions (/h^ 0.67, a;(CO) = 10"'', a;(0) = 3 x 10""^). At 
= 1.5 X lO"** equation (jl8p gives an ionization rate 1.33 times larger than equation (jl5p . so if anything, we are slightly underestimating 

C2. 
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Fig. 15. — Cosmic-ray ionization rates and upper limits from the sample of diffuse cloud 
sight hues. The left panel shows ionization rates inferred from detections of with la- 
uncertainties. Points are spread out horizontally for clarity. The right panel shows 3a upper 
limits to the ionization rate inferred for sight lines where is not detected. Some of the 
3(7 upper limits suggest ionization rates lower than those inferred in sight lines with 
detections. This likely points to variations in the cosmic-ray ionization rate between sight 
lines. 
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Fig. 16. — Upper panel: Probability density function, /(C2), found by combining inferred 
ionization rates and upper limits. This is defined such that the probability of the cosmic-ray 
ionization rate being in any particular range is given by P(a < (2 ^ b) — /(C2)f^C2- Center 
panel: Probability density function multiplied by the ionization rate, (2/ {(2)- This allows one 
to see "by eye" what portion of /(C2) carries the largest probability given the logarithmic 
X-axis. The vertical dashed line shows the mean value of (2 = 3.5 x 10~^^ s^^, and the 
shaded region is the smallest range of ionization rates in log space that contains 68.3% of the 
area under the probabihty density function. Lower panel: Cumulative distribution function, 
F{(2), showing the probability that C2 is below any given value. 
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Fig. 17. — A^(H3 ) versus A^(H2) for diffuse clouds studied herein and various dense clouds 
reported in the literat ure. Symbols a re as follows: filled diamonds - this study (detec- 
tions); filled squares - llndriolo et al.l (12010); gr ey arrows - this st udy (3a upper limits 



on A^(Hj")) ; open square s - 
triangles 



McCall et al. 



19991) 



open diamonds - iKulesal (|2002[ ): open 
Brittain et al.l ( l2004l ) and iGibb et al.l ( l2010l ). The diagonal lines show con- 
stant values of A^(H3 )/A^(H2) and are labeled accordingly. Diffuse clouds cluster about 
N{}l'^)/N{B.2) = 10~^, while dense clouds fall in the range between 10~^ and 10~^. For dif- 
fuse cloud sight lines A^(H2) is determined from UV H2 observations, estimated from A^(CH), 
or estimated from E{B — V), and values used here are presented in Table HI For dense cloud 
sight lines N{B.2) is determined from IR H2 observations, estimated from Ay, or estimated 
from N{CO). 
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Fig. 18. — Top panel: Map showing the locations of targeted sight lines in Galactic coordinates. Black diamonds 
mark sight lines where H;^ absorption is observed, and open grey triangles mark sight lines where Yl^ is not detected. 
Middle panel: Heliocentric distance of background sources with respect to Galactic longitude. Black lines mark sight 
hnes where absorption is observed, and grey Hnes mark sight lines where is not detected. Bottom panel: The 
cosmic-ray ionization rate as a function of Galactic longitude. Black diamonds are inferred cosmic-ray ionization rates 
with la uncertainties. Grey arrows are 3a upper limits. 




Fig. 19. — The cosmic-ray ionization rate as a function of cloud column density, N^, and equivalent color excess, 
E[B — V). For sight lines with multiple distinct velocity components, the total line-of-sight iVn has been divided by 
the number of cloud components. Filled diamonds are cosmic-ray ionization rates with la uncertainties inferred for 
diffuse molecular clouds in this study. Grey arrows are 3a upper limits from this study. Ionization rates inferred in 
dense cl ouds are as fol lows: open diamonds - observations of NGC 2024 IRS 2, AFGL 490, AFGL 2591, and AFGL 
2136 bv iKulesal fl2002h: crosses - HCO+ observations of AFGL 2136, AFGL 2591, AFGL 490, W33 A, and W3 IR S 
5 by Ivan der Tak fc van Dishoeckl (20001): open square - HCNH"'" observation of DR 21(011) by iHezareh et al.l ( 120081 ). 
Various models (e.g.. iPadovani et al.ll2009l ) predict that (2 should decrease with increased A^h as low-energy cosmic rays 
lose all of their energy in the outer regions of a cloud. While no strong correlation is apparent from the diffuse cloud 
data alone, the addition of the dense cloud ionization rates from sight lines with much higher A^h does suggest such a 
trend. 



